Steady-state solution for dark states using a three-level system in coupled quantum dots 
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Quantum dots (QDs) are one of the promising candidates of interconnection between electromagnetic field 
and electrons in solid-state devices. Dark states appear as a result of coherence between the electromagnetic 
fields and the discrete energy levels of the system. Here, we theoretically solve the steady-state solutions of the 
density matrix equations for a thee-level double QD system and investigate the condition of the appearance of a 
dark state. We also numerically show the appearance of the dark state by time-dependent current characteristics. 



I. INTRODUCTION 

Future quantum communication systems might be com- 
posed of optical fiber networks and quantum computers based 
on solid-state qubits. Even though quantum computers made 
by optical systems would be smoothly connected to optical 
communication networks, in order to interface to mobile elec- 
tronic systems, such as mobile phones, quantum computers 
based on solid-state circuits are desirablsiia. In this respect, 
an efficient interconnection between optical and solid-state 
systems should be developed. Quantum dot (QD) systems, 
such as GaAs/AlGaAs^^ have discrete energy levels which 
are suitable for the transfer of photon or phonon energies to 
electrons in solid-state circuits and can be used as elements of 
qubits'^'P. QD systems also have the advantage that the dis- 
tance between energy-levels can be controlled by the bias, in 
addition to the sizes of the QDs. Thus QD systems are one 
of the promising candidates for the interconnection between 
optics and solid-state circuits. 

One of the efficient connection methods is constructed by 
using coherent population trapping (CPT) or electromagnet- 
ically induced transparency (EIT)ii. CPT is a typical phe- 
nomena of quantum coherence in three-level system and has 
been intensively studied in opticsi^"— . By adjusting two laser 
fields, the electron population between the lowest two energy- 
levels are coherently transferred. Here, we theoretically dis- 
cuss the interaction between an electromagnetic field and a 
coupled double QD (DQD) system by focusing on transport 
properties of three-level systems. 

Recently CPT has been studied in superconducting qubit 
system&ii"— . Double quantum dot (DQD) systems such as 
GaAs/AlGaAs^— or Si/Si02 are also candidates for realiz- 
ing three-level systems and have been theoretically investi- 
gated22i^. Tokura et alJ^ theoretically investigated resonant 
tunneling currents under locally different Zeeman energies 
and found that when the magnetic fields in each QD are non- 
collinear, four resonant peaks can be observed. Ke et ai— 
constructed density matrix equations and investigates the re- 
lation between the phase of the driving lasers and transport 



properties. Emary et alJ^ investigated transport properties 
when three-energy levels exist in the same QD. From an ex- 
perimental viewpoint, applying two laser fields is not easy to 
control. Here, we mainly discuss the case where one of the 
laser fields can be replaced by electronic tunneling between 
two QDs. 

A three-energy-level DQD system is realized under a large 
bias voltage as depicted in Fig. 1, in which there is one energy- 
level (£1) in the left QD and two {E2 and £3) in the right QD. 
We assume a strong Coulomb interaction between electrons 
such that only one excess electron is allowed in the two QDs. 
We also assume that the left energy level Ei is close to the 
right upper energy level £3 such that electrons in QDl tunnel 
directly into £3 (£3 - £1 «: Q.^; and Q/, is the tunneling rate 
between QDl and QD2). 

The dark state is a state in which there is no electron in 
the £3 level and induces interesting phenomena in the trans- 
port properties of the DQD system. Let us first think about 
the case of conventional tunneling processes without the dark 
state: an electron tunnels from the left electrode to the QDl 
with a tunneling rate Fi. When the laser pulse is switched 
off, the electron is trapped at the £2 level with some proba- 
bility. Once the electron is trapped at the £2 level, because 
of the Coulomb blockade effect, there is no current through 
the DQD. When the laser pulse is switched on, the electron 
is excited from the £2 level to the £3 level. By the electron 
tunneling from the £3 level to the right electrode, the current 
finally flows through the DQD. Then, the Coulomb blockade 
is released, and a new electron can tunnel from the left elec- 
trode to the QDl. Thus, as long as the laser pulse sequence 
continues, the current continues to flow. However, when the 
dark state is realized, the three-level system is in a coherent 
superposition state, and because there is no electron in the £3 
level, the current does not flow through the DQD system in 
spite of the applying laser field and bias voltage. 

The purpose of this paper is to show the conditions for the 
production of a dark state, as functions of the Rabi frequency 
and the detuning parameter of the external laser field. We de- 
rive the density matrix equations in the three-level DQD sys- 
tem, and derive a steady-state solution for a dark state. We 
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FIG. 1: (a) A three-level system in double QDs (DQDs). A bias 
voltage is applied between the left and right electrodes, (b) A density 
matrix for the three level. We define = £3 - Ei and the detuning 
Ey = E^ - E2 - vjf. The dark state is a state with P33 = 0. 
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FIG. 2: The Rabi frequency Clu needed for realizing a dark state, as 
a function of the detuning Ey and yg. Here, = 0, 721 = 0, fo = 
and Fi = r2 = F3. These are solutions of P33 = in Eq.(|9}. 



investigate the relationship between the time-dependent cur- 
rent characteristic and the dark state. 

This article is organized as follows. In §2, the formulation 
of our model is presented. In §3, we show the analysis of the 
steady-state solution of the dark state. Section 4 is devoted 
to the numerical calculations for the time-dependent current 
when there is a dark state. The conclusions are given in §5. 
The detailed analytical solutions of the dark state in its steady 
state are shown in the Appendix. We also argue the analysis 
of the P22 - state in the Appendix. 



II. FORMULATION 



The Hamiltonian is H - Hq + H, + H; + Hy, where 
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£i|lXl|+ £212X21 + £3|3X3|, 
-(Qz.|lX3| + t^Re-''''"|2X3|) + h.c., 
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(1) 



Here |/) (/ = 1,2,3) is the energy-level in QDs, |/t/,> (\kR}) 
is the left (right) electrode state. is the Rabi frequency 
between E2 and £3, induced by an external laser field, and vr 
is the laser frequency. Vl (Vr) are the tunneling strengths of 



electrons between the left (right) electrode and the left (right) 
QD. The density matrix p,y = at T = is derived using 



dp i 1 

^ = -^[H,p]--{r,p), 



(2) 



where y expresses the dissipation of the system. In the present 
case of a QD system, the electron tunnelings between the QD 
and the electrodes correspond to this dissipation (hereafter we 
set /i = 1). When we define 
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the density matrix equations are given 



dt 

dpn 
dt 



= rf 'pH + Yfp22 + P33 - [Tf + r^^' + vf) 



Poo, 



^n^'poo - r'^'^pii + iQ.*J)3\ - i^Lpn 

+ito(pi2 -P2l), 

dp22 ,^(e)~ Wh)~ , n,* ~ ~ 
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i^lPii 



H-;QrP23 - in*RP32, 
= -(;w3i -I- y3i)P3i - i^L(P33 -Pn) + i^Rpii, 

= - ['(W32 - Vr) + 732] P32 

-/Qr(P33 -P22) + i^Lpn, 
= - [;(w2i + Vr) + 721] P21 - i^LP2i + i^*RP31 
-itoipu -P22), 



(7) 



where 731 , 732 and 721 represent decoherence, such as acoustic 
phonons. In the present conditions, p,o and po, are given by: 



dpi 
dt 



-;(w/o + r/o)p/o, 



(8) 



These equations are solved analytically, but they are indepen- 
dent of the density matrix equations Eqs. (6), therefore, irrel- 
evant to main transport properties; ojij = E, - Ej (i =1,2, 3); 
rj" represents an electron tunneling from the DQD to the elec- 

(h) 

trodes, and F) represents that from the electrodes to the 
DQD, where = 27rp,(£f,)|y|'*|2, with pa(Efi) (i = 1,2,3 
and X = e,h), for each electrode at the Fermi energy Epi 
(Ef\ - iiL,Ef2 = Efj, = /ir). 

Depending on the relative positions of £1, E2 and £3, we 
can classify the electron transport into the following two re- 
gions. 

(1) rf = 0(£2»/^^), 

(2) F*'-' = (£2 « t^R). 
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Because there is a finite bias between the left and right elec- 
trodes, the electron does not flow into £3 from the eight elec- 



(a) 



trodes, such that F 



and r 



■(e) 



are satisfied. 



In this 

■(f) 



paper we consider the case shown in Fig. 1 and set - 0' 
Fj'"' - 0, and F*^*^' - (this is the case in which a dark state 
explicitly exists). Hereafter we consider 731 - yo - Jyi- 



III. STEADY-STATE SOLUTIONS FOR THE DARK STATE 

Steady-state solutions are obtained from the density matrix 
equations when dpijjdt — 0. Compared with the optical three- 
level2a, the existence of the poo state complicates the equa- 
tions. A dark state corresponds to the case where there is no 
electron state in £3 as P33 - 0. For a given DQD, we can con- 
trol the electron tunneling by adjusting the laser field (Qr and 
vr). When fo - 0, we can express the steady-state solution by 
fourth-order polynomial equations of E^, Ey and Qr, such as 



P33 ^ A33(Or/Ol)^ + B^mRl^Lf + C33 = 0, 



where 



A33 = r;(-r^ + r2i). 

B33 = -E'yD'r^y',,+Dl\Y2' +yo')yo' 
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(11) 



C33 = r2'ro'[(l + D[E';)- + yl.'El' + (Of + yl^')yl' 
+ 2y'2,y',], (12) 

with D. = £1 - Ey = E2 - El + vr (all quantities are 
rescaled by and indicated by the prime symbol, such as 
F^ - T2/0.l)- This equation is a parabolic function regard- 
ing with C33 > 0. Thus, if F2 > 721 > P33(f^R) - has a 
solution for positive Q^. Also when F2 > 721, the coefficient 
of has a negative value, therefore, the fl^ of Eq.(|9]l for the 
dark state is a maximum value for the solutions of the density 
matrix equations to be valid. Figure |2] plots Qr, which satis- 
fies P33(0.l) = as a function of Ey/D.i. It can be seen that 
the larger Qr is required as Ey or decoherence yo increases. 

Equation ( fT2] i shows the relationship between and £2 - 
El + vr. Figure[3]show £,, as functions of £2 - £1 + vr and yo 
for 721 =0, and 721 = O.SQ/,- We can see that the region of the 
existence of £,, for the dark state becomes smaller in particular 
when 721 becomes larger. When 721 =0, £,. is given by 

£v = -^-^(-VrI+ ^F2Q^-D2(F2+7o)a^-Z)?-r27^) 

(13) 

Because the equation in the root square in this equation should 
have a real solution, we have the condition; 



£2 - £1 Vr > 



2F2yo 
2F2 + 70 



(14) 



Thus, when the decoherence 70 is larger than the tunneling 
rate F, £2 - £1 -n Vr > 2F, and when the decoherence 70 is 
smaller than the tunneling rate F, we have £2 - £1 + vr > 70. 
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FIG. 3: Condition of the dark state of Ey as functions of E2 — E1+ vr 
and 7o when (a) 721 = O.Sfit and (b) 721 = 0. Eq.([9]l is solved for Ey. 
We take to = 0, fl/; = 2f2z. and T2 = 



IV. TIME-DEPENDENT CURRENT 

Here we show numerical results for the time-dependent ma- 
trix element P33 and current. Qr is calculated from Eq. ^ 
such that the initial £,. is given, e.g., as £v/Ql = 2. FigureH] 
shows the time-dependent density matrix element P33 when 
£1 = 0, fo = 0, and 721 - 0, starting from (a) |1) and (b) 
(|1) + |2))/2. It can be seen that, as Ey decreases, p33(f) de- 
creases. 

As mentioned above, Qr is determined such that it satisfies 
the steady-state solution p33(f — > c») — > for £y/Q£ = 2, 
where p33(f — > 00) has the lowest values. Compared with 
Fig. lUa), Figure Hfb) oscillates faster. This is because, for 
the superposition state, the density population of electrons os- 
cillates between |1) and |2) more often than the case starting 
from |1). 

Figure |5] shows the time-dependent currents through the 
DQD system. The current is derived^^ as 



/(f) = FR[p33(f)+P22(f)]■ 



(15) 



Here we consider /(?) s e""'I{t). Qr is determined simi- 
larly to Fig |4] Thus the current is expected to be reduced for 
EylQ.1 - 2. Figures |5|a,c) show that the current decreases 
around the expected dark state. Figures [Sjc, d) show the 
time-dependent currents starting from a superposition state of 
(|1) + |2))/2. Compared with Figs.|5ja, c), Figs|5jb, d) show 
that a finite leak tunneling (fo - 0.5 and 721 - 0.5) leads to 
a small current reduction, and the evidence of the dark state 
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FIG. 4: Time-dependent density matrix elements pij,{t) starting from 
an initial state of (a) |1) and (b) (|1) + |2»/2. Here, yalSlt = 1, fo = 
and 721 = 0. Also, Ey/Cli = 5 corresponds to the solution of Eq.([9]l. 



disappears regardless of the initial state. 



I/e 



f/n, 
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V. CONCLUSIONS 

We theoretically solved the steady-state solutions of the 
density matrix equations for a thee-level DQD system, and 
showed the condition for the appearance of a dark state. Nu- 
merical calculations for time-dependent current characteris- 
tics showed that the steady-state can be detected by measuring 
a current. 
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Appendix A: Explicit form of density matrix 

Because the equations for fo = are still too complicated, 
we describe the case of - 0. 

pii = £{roi^i'((£v' + r2i')ro(i+r'ro) 

+ [£f + r2i(2r' + r^i)]r^'(i+Q2') 

+ 4'Q^'(r2i + r'ii^') + ro[<' + (-1 + + ^1') 
+ Ef(-2 + rl' + 2nl' + 2r'r2iQ^')] 

+ r^'[2r^i(i + Q^') + r'/i2]) 



P22 



(Al) 
(A2) 



FIG. 5: Time-dependent current as a function of starting from 
|1> for (a) and (b), (|1) + |2))/2 for (c) and (d). Ft/Qi, = Tj/n,. = 
Tj/ai. = 1, yo/fiz. = 1, (a)(c) to/Slt = and yji/fli. = 0. (b)(d) 
to/ilt = 1 and 721 /^L = 1- 



P33 = -|r(,Qpii+r^(2y^iro' + (£v' + r2i')ro' 

+ <r2i"«' + ro[i 

+ <(-2 + r?i' + n^')])|, 

Zo = r2'[ro'f^«'«(-i+ro') + (i+r2iV + "^)') 
+ r2'yo^iin+r^{Ef + yl')yt' 

+ E^:^l'; + 2y,,yl'{l+Q:l') + yl'h2) 



(A3) 



+ r 



4y,nl'h, + r',{3Efy,,iil' 
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r^iro'(6 + ^1' + 2Q^') + ro(3 + 3£f - 20^' 

Q^' + 3£f(-2 + r2/ + Q2')) 



where 



/12 = <'+<(-2 + r^/) + (l+Q2')2. 



(A4) 



(A5) 
(A6) 



Here, all quantities are scaled by such as = 0.^/0.1, 
r' - r/Qi and so on. 



+ [2ro(r-ro) + rr2i]£v£v 

+ yoirii + n [q^ + 72170 + n^] = 



Thus, Ey is written as function of 



where 



(Bl) 



Ey = ^ {Z72£.v + ^(/72 - 4fl2)£2 + 4^202} (B2) 



Appendix B: P22 = state 

In the main text, we discussed the dark state condition of 
P33 = 0. Here, we consider the region of p22 - without 
considering the dark state. Because there is no direct tunnel- 
ing term between the Ei level and the E2 level, an electron 
exists at the energy-level E2 only when there is some relax- 
ation process of the electron from the £3 level or dissipation. 
If the electron is transferred from the left electrode to the right 
electrode without the £2 level, we can regard the system as 
a two-level system as if each QD has one energy-level. The 
condition p22 - corresponds to: 



02 = 7o(r-7o) (B3) 
b2 = a2 + Ty2i (B4) 
C2 = 70(721 +r)[Qi + 72170 + (B5) 

This is the condition that the present system can be treated 
as a two-level system. When 721 - 0, we can simplify the 
condition as: 



Ey = £, + 



ml + ^l) 



\ r-70 



(B6) 



P22 ^ -7o(r-7o)(£', +£■;) 
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